Complete transfer of populations from a single state to a pre-selected superposition of 

states using Piecewise Adiabatic Passage. 
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We develop a method for executing robust and selective transfer of populations between a single 
level and pre-selected superpositions of energy eigenstates. Viewed in the frequency domain, our 
method amounts to executing a series of simultaneous adiabatic passages into each component of 
the target superposition state. Viewed in the time domain, the method works by accumulating the 
wavefunction of the target wave packet as it revisits the Franck Condon region, in what amounts 
to an extension of the Piecewise Adiabatic Passage technique [ Shapiro et.al., Phys. Rev. Lett. 
99, 033002 (2007)] to the multi-state regime. The viability of the method is verified by performing 
numerical tests for the Na2 molecule. 
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I. INTRODUCTION 

The availability of robust and selective methods of ex- 
ecuting population transfers in multilevel quantum sys- 
tems is essential for a variety of fields, such as preci- 
sion spectroscopy p , H, Q|, [Bl, quantum computing 
@, 0, H, B El, [lH, control of molecular dynamics 
and chemical reactions [l^, [lH , production of cold 
molecules [l6|, [13, [HI , biophotonics [l^ and nanoscience 
[2Q| , [21|. In this paper we propose a new method for 
transferring populations from a single energy eigenstate 
into a selected superposition of states (wave packet) us- 
ing shaped broadband laser pulses. This simple method 
combines the robustness of adiabatic population passage 
[22I . [23I . [24] with the flexibility of femtosecond pulse shap- 
ing techniques [13„ J^, ^5, 26]. 

Our method is an integration of a number of ear- 
lier studies. Viewed in the frequency domain, it is an 
application of Coherently Controlled Adiabatic Passage 
(CCAP) [21I, [53], which in itself is an extension of the 
three-state Stimulated Raman Adiabatic Passage (STI- 
RAP) [24, 28, 29] and the adiabatic transfer between 
field-dressed states 0, 0, 0, |^, 0, Q methods. The 
CCAP method presents a complete solution to the non- 
degenerate quantum control problem, i.e., the execution 
of a complete population transfer between superpositions 
(wave packets) of non-degenerate energy eigenstates. 

Viewed in the time domain, the present method is an 
extension of Piecewise Adiabatic Passage [13, [H, [13 to 
the case of a target wave packet of states. In particular, 
when the spectrum of the target manifold is nearly har- 
monic, the driving field is given by a train of mutually 
coherent pulses separated by the evolution period of the 
target wave packet. Keeping the optical carrier phase of 
the pulses constant throughout the entire train of pulses 
results in the piecewise execution of periodic (Rabi) os- 
cillations between the initial and the target states. The 
introduction of a "piecewise chirp" , expressed as a pulse- 



to-pulse variation in the optical phase, is what eliminates 
the oscillations and renders adiabatic robustness to the 
population transfer. 

Selectivity is obtained by tailoring the temporal and 
spectral profiles of the train of pulses to the attributes 
of the target wave packet dynamics |38|, |39|, |40| . Such 
compatibility between the pulse train attributes and the 
target wave packet dynamics has also been noted in 
brute-force optimization studies [is', 'ZT, 42] aimed at 
either maximizing population transfer into the target 
state [43I, [13, [in or stabilizing such transfer against wave 
packet spreading and decoherence [46^. 

The time profile of the field in our solution is reminis- 
cent of the "multi-RAP" pulse sequences of Ref. [35] . As 
discussed below, the difference between the two solutions 
is manifest when the target wave packet consists of more 
than two eigenstates. Our method is also related to the 
coherent accumulation of transition amplitudes driven by 
a train of laser pulses U, i, 113, IH, [HO, [III, [H , but 
the robustness is a unique property of PAP. 

This paper is organized as follows: In Section [TTl we 
present the theory in the frequency space and illustrate 
it via numerical studies of population transfer in sodium 
dimers. In Section HII} we describe the population trans- 
fer in the time domain, and connect the two physical 
points of view. In Section IIV| we establish numerically 
the accuracy of the method and its sensitivity to the pulse 
parameters. A discussion is provided in Section [Vl 



II. MULTI-STATE ADIABATIC CHIRPING - 
FREQUENCY DOMAIN POINT OF VIEW. 

In this section we extend the CCAP |2l|, [53] method 
of using a multi-mode pulse to execute wave packet adi- 
abatic passage via an intermediate state to the chirped 
multi-mode pulse case. As in Refs.0, consider a 

quantum system, initially in the ground state |0), coupled 
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by a laser field field e(t), made up of discrete modes, 
of frequencies a;^, to a manifold composed of N excited 
states. Contrary to Refs. [2l|, [l^ we allow the mode fre- 
quencies to slowly vary with time, hence we denote them 

as COnit). 

The material Hamiltonian is 

N 

n 7 

n=l 

where Hm is the field- free Hamiltonian, 2f{t)en and <pri 
are the mode amplitudes and phases, respectively, ft is 
the dipole moment for the transition between the ground 
state and the manifold of excited states. We assume that 
each mode frequency conit) = ^n{t) is detuned by a small 
amount An{t) = En — Eq — U0n{t) with respect to one of 
the excited levels (denoted n), where Eq and E^ are the 
field- free energies of the states |0) and |n), respectively. 



(^0 - HmM = {En - HM)\n) = . 



(2) 



The material wave function expressed in a.u. {h = 1) is 
expanded as 



N 



n=l 



It is instructive to study the dynamics arising from the 
application of the Rotating Wave Approximation (RWA) 
[li, H, [13, m. To do so we first average the non- 
stationary Schrodinger equation over time scales of the 
order of l/cjn, resulting in 



N 



N 
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(4) 



1 = 1 



N 



ibn{t) = A{t)bn{t) ^ f{t)bo{t)flnoYeii 



,i((uJn-UJl)t-cl)l) 



1 = 1 



We also assume that each transition is driven by 
only one of the field modes - the mode with the 
smallest An{t) detuning. Denoting by ftn{t) = 
/(t)en/inO e"*"^"" the (complex) Rabi frequencies, and by 

h{t) = {bo{t)^bi{t)...bn{t))^ , the (column) vector of ex- 
pansion coefficients, with the superscript T marking the 
"transpose" operation, we can now write the Schrodinger 
equation in matrix form as 



where 



zb = Hb 



/ ni{t) 



(5) 







AN{t) J 



(6) 



While the first part of the RWA (Eq.g])) amounts to 
neglecting the terms which oscillate at optical frequen- 
cies, the second part (Eqs. (|5|6p ) is equivalent to averaging 
the Schrodinger equation over time scales of the order of 
l/{u;n — ^m), resulting in loss of information about times 
scales shorter than a vibrational period. Nevertheless, 
the time-averaged Schrodinger equation ([5j) provides an 
accurate description of the wave function at the end of 
the process, after many vibrational periods. We post- 
pone the study of the dynamics on a finer time scale to 
the following sections. 

We now tune the chirping of the mode frequencies 
such that all the detunings are equal to one another, 
i.e., An{t) = A{t). This allows us to easily diagonalize 
H at any given moment of time, obtaining — 1 degen- 
erate "dark" eigenstates whose (quasi-) energies are equal 
to A(t), and two "bright" eigenstates whose eigenvalues 
are. 



X±{t) = A{t)/2±{A^{t)/4 + nls{t)P 



(7) 



where ^eff(^) ~ Xln=i l^n(^)P- The eigenvectors corre- 
sponding to the bright solutions are 



D± = cos u± 



sin 0±hf 



(8) 



where 



b, = (i,o,...,o)T bJ = (o,^^l,...,^],v)^Mff , (9) 

with tan^i = A±/r^eff- Note that each Qn^ as well as 
l^eff, depends on time via the common factor f{t). There- 
fore the definition of bj in Eq.([9]) is time independent. 
Further, tan^+ • tan^_ = — 1, and A_ < 0, hence. 
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If at some instant \A{t)\ ^ r^eff(^), then one of the states 
represented by the b+ or b_ vector coincides with the 
I i ) state, represented by the b^ vector, and the other - 
by the | /) state, represented by the b/ vector. 

Similar to the textbook case of adiabatic following in a 
two-level system [23], when A{t) is made to vary slowly 
enough, the system makes a smooth transition from state 
I z ) to state I / ) . We show in the Appendix that there are 
no transitions between the | ± ) bright eigenstates and the 
dark states. As a result, the dark states remain unpop- 
ulated at all times. Hence it is possible to make the adi- 
abatic population transfer complete, provided that the 
bandwidth of the pulse covers the whole range of the 
target energies En. 

In order to create the required field, one can start with 
a single broadband laser pulse, and spectrally shape it 
in the way illustrated in Fig. [H^a). Here, the field is 
frequency chirped in the neighborhood of each resonance 
frequency co^^^ = {En — Eq)^ according to. 



gi(wto-<^>„) 



(11) 



The real amplitude envelope F{lj) reaches its maximum 
of 1 at = and serves to suppress the pulse spectrum 
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FIG. 1: (Color online) Piecewise adiabatic following in 
Na2. (a): Amplitude spectrum (solid black) and frequency- 
dependent phase (dashed blue) of the driving field, (b): The 
driving field as a function of time. (c,d): Populations of var- 
ious vibrational levels of the ground and excited electronic 
states, (e): Final populations of A^T^i states predicted by 
the theory (full red circles), and obtained in the numerical 
simulation (empty blue squares). (f): Same as (e) for the case 
when the spectral region corresponding to the transition into 
one of the excited levels is blocked. 



between the resonance regions. The requirements for this 
suppression will be discussed in Section lTVl For Gaussian 
profile of the field envelope near each cc;J^^*, 



F{u, - cv^n = exp [{iv - <^)V2<tS] . (12) 



The above form corresponds to the field parameters of 



Eq. dll) assuming the form, 

m 



o-™(27r)i 



at 



-arg [1 - ia^al] /2 



rexp[-{t-tof/2a^] 



(13) 



In the adiabatic transfer into a pre-selected superposi- 
tion state, the real amplitudes = e(cc;J^^^) and phases 
(j)n are chosen such that the vector b / given by Equation 
([9]) represents the | / ) target state. The direct corre- 
spondence between Eqs. (|7|8p and the equations describ- 
ing adiabatic following in a two-level system driven by 
a single-component chirped pulse [23*] suggests that the 
variation of each chirped resonant frequency must exceed 

Fig. [T] shows a simulation of the adiabatic transfer be- 
tween the = vibrational state of the Na2(X^S+) 
ground electronic state, and a vibrational wave packet 
composed of the v = 7, ..,12 states of the A^S^ state 
[53^ . In order to assess the fidelity of the transfer we have 
also computed the population of neighboring X^T^'^ and 
vibrational levels. Panel (a) shows a field spec- 
trum obtained by frequency modulating a sin^ pt-type 
pulse. The pulse duration is 55fs and its central wave- 
length is 638 nm, with chirp parameter a^j = 2 x lO^fs^. 
As discussed in detail in the next section, the piecewise 
chirped pulse of panel (a) corresponds in the time domain 
to the train of pulses shown in Fig. [TJb). 

In Fig. [ifc-e) we display the dynamics of the popu- 
lation transfer. As demonstrated in Fig. [TJe), the pop- 
ulations of the target wave packet energy components 
practically coincide with the predictions of Eq.([9|). In 
Fig. [TJf) we demonstrate the selectivity of the method 
by changing the amplitudes and the phases (j)n char- 
acterizing the field components. In each case, the dis- 
tribution of the final populations in the target manifold 
remains close to that predicted by Eq.(|9|), though the 
component phases could at times deviate from the target 
ones, somewhat lowering the magnitude of the overlap 
with the target wave packet to ~ 0.8 — 0.95. This point 
is discussed further in section [IVl 

The mechanism of the population transfer becomes 
transparent if the laser field is represented by a phase- 
space (Husimi) spectrogram. Fig. [2fa) shows the spec- 
trogram of the field of Fig. [Ha,b). Each point displays 
the absolute value of the overlap integral of the field with 
a Gaussian probe pulse. 



e^^^(t) =exp 




(14) 
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of long duration {a^^ = 1000 fs). The frequency and 

time axes correspond to u^^"^ and respectively. It is 
clear that the field drives a number of adiabatic passages 
simultaneously: An adiabatic transition into n*^ excited 
level may occur when the instantaneous frequency com- 
ponent ujn crosses the n*^ resonance at uj^^^^ = En — Eq. 

The mechanism displayed here is quite different from 
that of the "molecular 7r-pulse" [13, [13, [H, [13, [58|] , where 
a transition into the manifold of excited states is driven 
by a single frequency-chirped laser field. While in the 
scheme proposed here all the transitions take place si- 
multaneously, in the case of a molecular 7r-pulse the en- 
ergy of the initial state shifted by the energy of one pho- 
ton, [£^0 + ^{t)]^ crosses all the target energies one by 
one. As a result, the system passes through a chain of 
avoided crossings. This is illustrated in Fig. [2jb) which 
shows the spectrogram (a^.^^ = 200fs) of the molecular tt- 
pulse obtained from the original un-shaped field, used in 
our numerical examples, by applying a single frequency 
chirp of = 2 X 10^ fs^. Although a molecular 7r-pulse 
can transfer all the initial population into the excited 
manifold in a robust way, the relative amplitudes of the 
states in the resulting wave packet are not controllable 
unless the target manifold consists of only two states. In 
the latter case, employing either up- or down- frequency 
chirp allows population transfer in either the lower or 
the higher eigenstate of the target manifold [s^. We re- 
fer the reader to Ref.[58] for a detailed analysis of the 
population transfer by molecular 7r-pulses. 

The question of selectivity also arises when compar- 
ing our method to the family of Stark-assisted adiabatic 
transfers [59], and the transfer via multiple successive 
Rapid Adiabatic Passages [35]. While in both examples 
the possibility of robust and selective transfer into either 
one of the two target states has been found, it is not clear 
whether these methods can enable population transfer 
into a selected superposition of target states, especially 
if AT > 2. 




Time (ps) ^ ° Time (ps) "'' Time (ps) ^ 



FIG. 2: (Color online) Field spectrograms, (a): Absolute 
values of the overlap integral of the pulse train field with 
probe Gaussian pulses. The red horizontal lines show the res- 
onance frequencies a;^^*. (b): Spectrogram of a conventional 
frequency-chirped pulse, (c): Expansion of the pulse train 
field in short probe pulses. The amplitude of each projection 
is shown by brightness, its phase by color. The white dotted 
line marks the spectral center of the pulse. 



III. MULTI-STATE ADIABATIC CHIRPING - 
TEMPORAL POINT OF VIEW. 

In order to understand the dynamics of the population 
transfer on a finer time scale, here we develop the time- 
domain description of the process. Fig. [H^b) shows the 
time-dependent driving field of the above example. The 
time analogue of Eq. (pT]) is 

e{t) = Re f{t) exp [iuo{t - to) + iat{t - tof/2 + ^i(^) 

(15) 

where ujq = cjj^^^ is chosen to coincide with the frequency 
of a transition into one of the central states |no) of the 
excited wave packet, and 

N 

^i(t) = 2^enexp[z(^n-^nJ(t-to)+i(/>n] . (16) 

n=l 

Equation (p!6|) is similar to an expression describing wave 
function dynamics of a fictitious wave packet in the tar- 
get manifold, with the amplitudes of eigenstates given 
by enex.p[i(j)n]. If the target spectrum and the distribu- 
tion of the complex field amplitudes exp[z(/)n] are both 
smooth, then the interference between the different fre- 
quency components of the field results in a train of short 
pulses separated by the vibrational period TyUj = dEn/dn 
of the target wave packet. The field pattern is restored 
whenever {En — Eno)t :^ 27rm for all n and any integer m. 
Such a train of short pulses is seen in Fig.flfb). Further, 
if the target spectrum is harmonic, then the field spec- 
trum is a frequency comb with equally spaced teeth, and 
the pulse shape is preserved within the train. If the spec- 
trum is anharmonic, then the shape evolves from pulse 
to pulse within the train, refiecting the spreading of the 
target wave packet as it revisits the transition region. 

When the target spectrum is weakly anharmonic, the 
pulse train of Eq. (p!5|) is composed of pulses which dif- 
fer from one another in four essential ways: First, the 
temporal spacings between the adjacent pulses (which 
may have different overall power) change. Second, their 
central frequencies differ by at Tyi^. Third, and most im- 
portant, the l^^ pulse has an extra overall phase equal to 
at {ITyii^Y where / = corresponds to the middle of 
the pulse train; and fourth, owing to the anharmonicity, 
the condition {En — Eno)t = 27rm is fulfilled for differ- 
ent n at slightly different times, leading to change in the 
pulse duration from pulse to pulse. 

Fig. [Sfc) shows the spectrogram of the field of 

Fig. [Ha,b) using short (a^^^^ = 50fs) probe pulses, whose 
exact form is given by Eq. ([Hj). Naturally, the repre- 
sentation in terms of short probe pulses gains temporal 
resolution but loses frequency resolution. The optical 
phase of the short pulses within the train varies quadrat- 
ically with the pulse number. This is displayed in Fig. 
[2](c) by the variation in color along the central-frequency 
section, from orange (phase equal to zero) in the mid- 
dle of the train, to purple and blue, and finally to green 
(phase :^ 37r/2), at the very tails. Note that the pulse- 
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to-pulse drift of the central frequency is not noticeable 
on the scale of the figure. 

The field in the above example is similar to that in- 
troduced in Ref.[37] for the piecewise adiabatic passage 
in a two-level system. This is readily seen both in the 
frequency-domain representation (by comparing FigUJa) 
with Fig. 3(b) of Ref.[37]), and in the time-domain rep- 
resentation. The almost-periodic pulse train only drives 
transitions between the | i ) and | / ) states. The train of 
short pulses, tailored to fit the periodic evolution of the 
I / ) state, peaks every time this state revisits the Franck- 
Condon region. At these instants a superposition state 
I ip ) can be represented by the unit Bloch vector of Fig.[3l 
defined by the and (j) angles, 

\i^)= cos{e/2)\i) + sin(^/2)e'*|/) , (17) 

The present 1+many- level system differs from the two 
level system in a fundamental way, since here the rela- 
tive phases of the components of the target wave packet 
change all the time. There are points of similarity as 
well, because here too, the process can be depicted, as 
shown in Fig. 3, by the discontinuous motion of a unit 
vector on a ("Bloch") sphere, with the polar angle ^ = 
representing 1^) and 6 = tt representing |/). The ef- 
fect of each pulse (of short duration r) is now viewed 
as a rotation P = R{ap) of the Bloch vector by an 
angle ap = J^Qefi{t)dt about the y axis. The change 
in the carrier phase between consecutive pulses can be 
represented by an additional rotation about the z axis, 
F = Rz{aF)' Thus, the overall evolution due to the 
pulse train of Eq. (p!5|) is represented by a sequence 
U = ...TFPTFP...^. The operator T, which corresponds 
to the free evolution of the wave packet between pulses, 
keeps the Bloch vector orientation at the beginning of 
each short pulse in the train equal to that at the end of 
the previous one. 

The product FP of two rotations can be viewed as an 
overall rotation by an angle ao about an instantaneous 
axis defined by the {Oo^(j)o) angles, given to lowest-order 
expansion in ap, as 

ao = {{ap^ + o^F^)/2}^ (^o = ±7r/2 - aF/2, 

tan6>o = iap/ap. (18) 

By maintaining the same value of ap and ap throughout 
the pulse train we can induce piecewise rotations of the 
Bloch vector about the "adiabatic" trajectory traced out 
by (^o,0o)- As shown in Fig. 3, by slowly varying the 
values of ap^ap we can make the Bloch vector follow 
this adiabatic trajectory. Such piecewise following can 
be realized provided, (i) the y- and z-rotations are small 
(i.e. each pulse should induce an angular change much 
smaller than tt and each increment in the carrier phase 
should be small too); (ii) (6>o, 0o) should not move much 
from pulse to pulse, i.e. 

A6>o<C{(ap2+ap2)/2}i (19) 



In the pulse train of Eq. ([15]), initially ap <^ |ap|: r^eff 
is small, while the pulse-to-pulse phase change is signif- 
icant. As ap increases and |ap| decreases, the states 
originating in \i) and |/) move towards the equator of 
the Bloch sphere. They cross the equator as soon as ap 
changes sign, and finally interchange with each other. 
Depicted in the original non-rotated frame, the trajec- 
tory of the Bloch vector is a piecewise spiral, similar to 
that shown in Fig|3l 



IV. NUMERICAL STUDIES. 

Fig. m shows the population transferred into the tar- 
get manifold, and the projection of the final state onto 
the target wave packet comprised of vibrational eigen- 
states of Na2 in the electronic state as a function 

of the chirp a^ and the field strength of the original 
55fs pulse used in the computations presented in Fig. [H 
The plots reveal several interesting features. First, for 
both positive and negative chirps exceeding in magni- 
tude |a^| ~ 200000fs^, the transfer is almost complete 
and quite robust with respect to changing a^j and e^. 
For a large range of pulse parameters the transfer proba- 
bility is ~ 95%, and the projection of the final state onto 
the target is - 0.85 - 0.9. 

The transfer probability landscape at small values of 
corresponds to the piecewise Rabi oscillations between 
the initial state and the excited wave packet. A close 
look at the population dynamics in this region of param- 
eter space shows that the first few pulses in the train 
manage to completely empty state |0), while the follow- 
ing pulses re-populate it. At stronger fields, state |0) is 
de-populated and re-populated several times during the 
pulse sequence. Populations of the nearby vibrational 
eigenstates of the X^T^^ manifold remain negligible at all 
times, with the ratios of populations of different 
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FIG. 3: (Color online). A calculated sample trajectory of 
the Bloch vector (thick gray arrow) during the a piecewise 
AP process implemented with a train of 20 ultra-short pulses 
|33]. 
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FIG. 4: (Color online) Population transferred into the target 
manifold (a), and the projection on the target wave packet 
(b) as a function of the chirp auj and field strength 



states remaining close to the target values. 

We also note that though the parameter-space region 
where robust population transfer occurs includes pulse 
areas of ~ tt, the optimal pulses have areas larger than 
TT. The pulse area cannot be increased beyond a certain 
value because when individual pulses within the pulse 
train are able to drive a significant population transfer, 
condition (i) above breaks down, and the overall fidelity 
decreases. 

The fidelity profile of Fig. [H is quite robust with re- 
spect to changes in the spectral width in the vicinity 
of each resonance frequency oo^^^ . However, when cr^ 
approaches half of the distance between the neighboring 
resonance frequencies, the transfer fidelity exhibits rip- 
ples as a function of a^, eventually becoming unstable. 
Indeed, for the Gaussian amplitude modulation of the 
spectrum fEq. (|12ill"3|) ). and assuming that the chirp is 
strong so that aa;Cr^ ^ 1, one has 



CTt 



at l/a^. 



(20) 



and the instantaneous frequency of the n}^ component of 
the field at the end of the pulse train becomes 



^n{t = ^0 + 2crt) 



2atat 



(21) 



If cr^ ^ (^n+i ~^n^^)/2 then the instantaneous frequency 
of the n}^ field component at the end of the pulse train 
coincides with the adjacent resonance. In this case the 
dynamics can not be viewed as a set of independent par- 
allel adiabatic passages into different target states. 



V. SUMMARY AND DISCUSSION. 

In this paper we have developed a method for execut- 
ing robust and selective transfer of population from a 
single energy eigenstate to a pre-selected superposition 
of energy eigenstates. Viewed in the frequency domain, 
the method constitutes simultaneous transfer of popula- 
tion to all the energy eigenstates which make up the su- 
perposition state by a set of parallel adiabatic passages. 
Viewed from the time domain, the method amounts to 
using a train of pulses to accumulate wavelets which make 



up the target wave packet as it revisits the Franck Con- 
don region [H, [sl, [s^l • 

We have tested the method numerically by simulating 
transitions between a single vibrational eigenstate and 
a superposition of vibrational energy eigenstates of Na2 
in the state. The simulations confirmed the high 

selectivity and robustness of the method. 

Topics to be investigated further include: An accu- 
rate description of the transfer dynamics for a general 
anharmonic spectrum and in particular the description 
of the short-time dynamics of the transfer. This may be 
achieved using Floquet states [311,|32[ dressed by a multi- 
mode driving field. Also, although the present theory 
is quite accurate in predicting the ratio of populations 
transferred into the target wave packet, the reasons for 
the phase errors, resulting in less than perfect overlap 
with the target state, need to be investigated. 
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VI. APPENDIX: PROOF THAT THE DARK 
STATES REMAIN UNPOPULATED. 

This proof follows Ref.[2l|, where similar arguments 
were made in the context of CCAP. These arguments 
were based on the standard treatment of adiabatic trans- 
fer 113|,|22|,|23|]. 

It follows from Eq.(j6|) that each dark state is de- 
scribed at a given time by the vector of amplitudes 
^dark = {^^ bi±^ ...^ which is orthogonal to both 

= (1,0,. ..0)"^ and hf = (0,l]i/l]eff,...,^^ArMff)"^. 
Since the mutual ratio of different Rabi frequencies 
Qri does not change in time, we can choose for the 
dark states an N — 1-dimensional time-independent basis 
(0, e^^\ e^^)''", n = l..A^ — 1. Introducing the vector 
a of amplitudes of instantaneous eigenstates. 



h{t) = U(t)a(t) 



(22) 



where the columns of U are given by the normalized 
eigenvectors of H, and using the relation U'^'U = 1, we 
have that 



ia = U^HUa-iU+Ua . 



(23) 



The first term on the right-hand side of Eq. (|23|) governs 
the adiabatic evolution. The second term governs non- 
adiabatic transitions between the instantaneous eigen- 
states. Let us choose the order of adiabatic eigenvalues 
such that the adiabatic Hamiltonian U^HU has on its 
main diagonal (A+, A_, A, ...A), with N — 1 terms equal 
to A, and all the off-diagonal elements equal to zero. 
Then, using Eq.([TO|), we obtain 



7 



COS ( 



sm ( 



sin ^+ r^i/r^eff — cos^+ r^i/r^eff e 





(1) 



sin ^+ l^TV /^eff — COS ^+ l^TV /^eff 



\ 



(24) 



As a result of the above structure of U, the non-adiabatic 
matrix is given as, 

Utu = I -^0 " I . (25) 




We see that this matrix does not couple the bright 



and dark states. Hence the TV 
unpopulated at all times. 



1 dark states remain 
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